Available online at www.sciencedirect.com 


ScienceDirect 


JOURNAL OF 


www.elsevier.com /locate /jpowsour 


Journal of Power Sources 165 (2007) 244-252 


Exergoeconomic analysis of a vehicular PEM fuel cell system 
S.O. Mert*®, I. Dincer**, Z. Ozcelik?” 


è Faculty of Engineering and Applied Sciences, University of Ontario Institute of Technology, Oshawa, Ontario, Canada 
> Department of Chemical Engineering, Ege University, Izmir, Turkey 


Received 17 June 2006; received in revised form 28 November 2006; accepted 1 December 2006 
Available online 16 January 2007 


Abstract 


In this study, we deal with the exergoeconomic analysis of a proton exchange membrane (PEM) fuel cell power system for transportation 
applications. The PEM fuel cell performance model, that is the polarization curve, is previously developed by one of the authors by using the 
some derived and developed equations in literature. The exergoeconomic analysis includes the PEM fuel cell stack and system components as 
compressor, humidifiers, pressure regulator and the cooling system. A parametric study is also conducted to investigate the system performance 
and cost behaviour of the components, depending on the operating temperature, operating pressure, membrane thickness, anode stoichiometry and 
cathode stoichiometry. For the system performance, energy and exergy efficiencies and power output are investigated in detail. It is found that with 
an increase of temperature and pressure and a decrease of membrane thickness the system efficiency increases which leads to a decrease in the 
overall production cost. The minimization of the production costs is very crucial in commercialization of the fuel cells in transportation sector. 
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1. Introduction 


During the past decade, alternative fuels and power sources 
have received an increasing attention due to the increasing 
environmental pollution and the decrease of the fossil fuel 
resources. So a movement towards environmentally friendlier 
and more efficient power production sources both for stationary 
and mobile applications is of paramount importance and brings 
the fuel cells to the forefront. 

Fuel cells are electrochemical devices that convert chemical 
energy of the reactants directly into electric and heat with high 
efficiency. Also, high operating temperature is not necessary for 
achieving high efficiency as the fossil fuels since electrochemi- 
cal processes in fuel cells are not governed by Carnot’s law. High 
efficiency makes fuel cells an attractive option for a wide range 
of applications, including transportation applications, domes- 
tic electricity, heat production and even portable and mobile 
systems. 
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PEM fuel cells are considered as the one of the most promis- 
ing alternative power sources especially for sub-megawatt scale 
applications like light-duty transportation and considered as a 
potential replacement with the classical conventional internal 
combustion engines. A PEM fuel cell powered vehicle has a 
number of important advantages as high efficiency, quick start 
up, low operating temperatures, high current density and zero 
pollution. On the other hand, having high capital costs, prob- 
lems of hydrogen storage and lack of infrastructure is the main 
problem for the commercialization of the PEM fuel cells in trans- 
portation applications. In these conditions, any increase in the 
system efficiency will help accelerate the commercialization of 
fuel celled vehicles [1]. 

Exergy analysis appears to be a potential for system design, 
analysis and process evaluation and improvement. Whereas 
energy analysis is based on the first law of thermodynamics, 
exergy analysis is based on both the first and the second laws 
of thermodynamics. Both analyses utilize also the material bal- 
ance for the considered system. Analysis and optimization of 
any physical or chemical process, using the energy and exergy 
concepts, can provide the two different views of the considered 
process [2]. 

There are only few studies on the exergy analysis of PEM fuel 
cells for transportation applications. Hussain et al. [1] investi- 
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Nomenclature 


a 
Acell 
AC, 
C 
Ck 


membrane activity 

cell area (cm?) 

annual capital cost (US$ year” !) 
exergetic cost (US$ kW7!) 
capital cost (US$) 


CRẸF(i, n) capital recovery factor 


PW 


exergy (kW kg7!) 

exergy, (kW) 

Faradays constant (C) 

flow rate (kmol s7!) 

enthalpy (J mole™!) 

higher heating value (kJ mole~!) 
current density (A cm); interest rate 
exchange current density (A cm7?) 
irreversibility (kW) 

life time, year 

number of fuel cells in stack 
pressure (atm) 

present worth (US$) 


PWF(i, n) present worth factor 


resistance (Q) 

universal gas constant (J (mole Ky 4) 
entropy (J (K mole)~!) 

salvage value (US$) 

membrane thickness (cm) 
temperature (K) 

cell potential (V) 

power (W) 

mole fraction 

capital Cost flow (US$ s7!) 


Greek letters 


a transfer coefficients 

n efficiency 

Amem membrane water content 

u chemical potential (J mole!) 

Omem membrane conductivity (1 (Q cm)—!) 
E stoichiometric ratio, ratio of excess air and fuel 
Subscripts 

A anode 

act activation 

C cathode 

ch chemical 

con concentration 

FC fuel cell 

irrev irreversible 

in inlet 

mem membrane 

ohm ohmic 

out outlet 

phy physical 

rev reversible 


gated thermodynamic analysis on fuel cell power systems. They 
studied the effects of operating temperature, operating pressure 
and air stoichiometry on the system performance. In addition 
to this study, Cownden et al. [3] performed exergy analysis of 
hydrogen fuel cell power system for bus transportation. Both 
of these two works showed the importance of the thermody- 
namic analysis in determining the irreversibilities within the 
system components and how the operating parameters affect 
the system performance. Ay et al. [4] derived a model for the 
thermodynamic analysis of a PEM fuel cell and found that ther- 
modynamic irreversibilities in the fuel cell increase with a rise 
of membrane thickness and with a decrease of cell temper- 
ature and pressure. For a PEM fuel cell using methanol, an 
exergy analysis has been studied by Ishihara et al. [5]. It is 
found that the energy efficiency approaches the unity as the 
recovery rate of the waste heat from the cell approaches the 
unity. The exergy efficiency is found to be about 0.45 with a 
fuel cell operating temperature of 80°C. It was also found that 
the cell voltage should exceed 0.82 V in order to obtain the 
exergy efficiency of 0.5 or higher. Moreover, Kazim [6] con- 
ducted exergy analysis of a PEM fuel cell at specified operating 
voltages of 0.5 and 0.6 V. In the study, the exergy efficiency 
of the PEM fuel cell is investigated depending on the oper- 
ating temperature (T/To) and operating pressure (P/Po) at the 
ratios ranging from | to 1.25 and from | to 3, respectively. Ay 
et al. [7] have proposed an exergetic performance analysis for 
a PEM fuel cell to investigate the effects of operating temper- 
ature and pressure to the system efficiency and irreversibilities. 
It is concluded in the paper that exergy efficiency of PEM fuel 
cell decreases with a rise in membrane thickness and current 
density, and increases with a rise of cell operating pressure and 
with a decrease of current density for the same membrane thick- 
ness. Kazim [8] also studied a brief exergoeconomic analysis 
of a PEM fuel cell at various operating temperatures, pres- 
sures and air stoichiometries. This study extends his previous 
study to cover the investigation of the air stoichiometry rang- 
ing from 2 to 4 from the exergy perspective. He indicated that 
the exergy cost of the fuel cell can be improved by adopting 
any or a combination of higher operating pressure, inlet air 
stoichiometry or cell voltage that demonstrates a significant 
improvement in the exergy cost. Barbir and Gomez [9] showed 
that there is a strong relationship between the efficiency and 
economics of PEM fuel cells. Typically, the lowest efficiency 
is achieved at maximum power output. The results indicate 
that in the best case scenario the fuel cells can be produced 
at US$ 100 kW7! with 50% efficiency, and generate electricity 
at cheaper than US$ 0.08 kW~! if hydrogen can be supplied at 
US$ 10GJ~!. 

Here, the present study differs from the above said studies 
since this investigates the effect of the operating parameters and 
some design parameters in a broad range including the system 
components other than the fuel cell stack, such as compressor, 
heat exchanger, humidifiers, pressure regulator and the cool- 
ing system on the system performance. Also, the present study 
investigates the exergetic cost flows of the system equipments 
and the total cost for the power production as an exergoeconomic 
analysis. 
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2. Exergoeconomic analysis 
2.1. Exergetic aspects 


Exergy analysis is an effective thermodynamic method of 
using the conservation of mass and conservation of energy prin- 
ciples together with the second law of thermodynamics for the 
design and analysis of thermal systems, and is an efficient tech- 
nique for revealing whether or not and by how much it is possible 
to design more efficient power systems by reducing the ineffi- 
ciencies [2]. 

The exergy concept was introduced to overcome limitations 
of the energy analysis. The exergy expresses the practical value 
of any substance (or any field matter, e.g., a heat radiation), and 
is defined as a maximum ability of this substance to perform 
work relative to human environment. Also, as Rosen and Dincer 
[2] pointed out, energy analysis embodied first law of thermo- 
dynamics only identifies external energy waste and losses, on 
the other hand, exergy takes entropy into account too by includ- 
ing irreversibilities. As a matter of fact, exergy is defined as the 
maximum theoretical work obtainable as the system interacts 
with its surroundings and comes to equilibrium. Once a system 
is in equilibrium with its surroundings, it is not possible to use 
the energy within the system to produce work. At this point, the 
exergy of the system has been completely destroyed. The state 
in which the system is in equilibrium with its surroundings is 
known as the dead state. 

In order to calculate the exergy of a system, we must specify 
both the system and the surroundings. The reference environ- 
ment is used to standardize the quantification of exergy. The 
reference environment or simply the environment is assumed 
to be a large, simple compressible system. The reference tem- 
perature is assumed to be uniform at Tọ, and the pressure is 
assumed to be uniform at Pg. Also, it is assumed that the inten- 
sive properties of the environment are not significantly changed 
by any process. Therefore, the environment is modeled as a ther- 
mal reservoir at To. In the present study, the restricted reference 
environment is chosen as 298 K and 1 atm and the unrestricted 
reference state composed of the N2, O2, H2O and CO2 and Ar 
in the molar fractions 0.775, 0.206, 0.018, 0.0003 and 0.0007, 
respectively. 

In general, there are several types of exergy as physical, 
kinetic, chemical, potential, nuclear, magnetic, electrical, etc. 
Here, we consider physical (Epny), chemical (Een), kinetic (Ex) 
and potential (Ep) exergies only as 


È = Ex + Ep + Epny + Eon = mex + ep + ephy + ech) (1D) 


The physical exergy is equal to the maximum amount of work 
obtainable when the stream of substance is brought from its 
initial state to the environmental state defined by Po and To 
by physical process involving only thermal interaction with the 
environment. 


phy = (h — ho) — Tols — so) (2) 


where subscript 0 stands for the reference environment 
(restricted). 


The chemical exergy is the exergy component associated with 
the departure of the chemical composition of a system from that 
of environment. The chemical exergy can be defined in a molar 
basis as: 


ech = 5x j(Mi; o — 47,00) (3) 
J 


where x; is the mole fraction of the species j in the flow; j4;,9 the 
chemical potential of species j in flow evaluated at Ty and Po; 
/4j,00 is the chemical potential of species j in the flow evaluated 
in the reference environment (unrestricted) [1]. 


2.2. Exergoeconomic aspects 


Exergoeconomic analysis is a method combining both exergy 
analysis and cost accounting. The method provides a technique 
to evaluate the costs of inefficiencies and/or the costs of individ- 
ual process streams, including intermediate and final products 
[10]. 

Exergoeconomics is nowadays a powerful tool to study and 
optimize a power system. The application field is the evalua- 
tion of utility costs as products or supplies of production plants, 
the energy costs between process operations or of an energy 
converter. Those costs are applicable in feasibility studies, in 
investment decisions, on comparing alternative techniques and 
operating conditions, in a cost-effective section of equipment 
during an installation, an exchange or expansion of an energy 
system [11]. 

The cost balance of a system may be written as [12]: 


X (Èin,iCin,i) + Zequipment = 5 (out, Coui) + WCw (4) 


where Finj, Eout,i Cini and Cout, are the exergies and exergy 
costs of the streams entering and leaving the control volume, 
respectively. Žequipment is the annualized cost of the equipment 
inside the control volume. Cw is the cost of the work or the power 
of the equipment. Using Eq. (4) the costs for the each component 
of the system is determined. Also, the cost balance is applied to 
the overall system to calculate the cost of the produced power 
by the fuel cell system. 

In order to calculate Žequipment» the annualized (or levelized) 
cost method is used [13]. The algorithm of this method is com- 
posed of four steps as outlined by Kwak et al. [14]. The first the 
present worth (PW) of the investigated system is calculated by 
substituting the effect of salvage value, Sg „n (Eq. (5)). In the cal- 
culations, the salvage values are taken as 10% of the capital cost. 


PWk = Ck — Sk.nPWF(i, n) (5) 


With the help of the CRF (capital recovery factor) that is, a 
function of the lifetime ( years) and interest rate (i), the annual 
capital cost is found, which is used for the calculation of the 
capital cost flow of the present system (Żequipment): 


AC, = PW x CRF(i, n) (6) 


, bx AC 
X= a eee (7) 
3600 sh~~ x 8000h yr 
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Fig. 1. Flow chart of Ballards Xcellsis™ HY-80 fuel cell engine [28]. 


i(1 + i)” 
CRF = ———_——_ 8 
(spi= 1 (8) 
where PWF(i, n) is the present worth factor that is taken as 
7.0236 [15], @ is the factor for the operating and maintenance 
costs which is taken as 1.06 of the AC. The lifetime of the 
system is taken as 10 years. 


2.3. The fuel cell engine system studied 


The Ballard’s Xcellsis™ HY-80 Fuel Cell Engine is taken as 
an example system for model calculations (Fig. 1). This engine 
is used by many automobile producers like Ford in Focus Fuel 
Cell Vehicle (FCV) and MAN Fuel-cell bus [16]. The engine 
is lightweight, 68 kW hydrogen-fuelled fuel cell engine that 
offers automotive manufacturers the opportunity to develop their 
own zero-emission fuel cell vehicles. The Xcellsis™ HY-80 
fits beneath the floor of the vehicle without reducing the size 
of the passenger compartment. It simplifies vehicle integration, 
assembly and service. The hydrogen stored in tank and fed to 
the system after a pressure regulation depending on the system 
pressure. Then, the inlet hydrogen is mixed with the recycled 
unreacted hydrogen from the fuel cell output. The final hydrogen 
stream is humidified in order to achieve the water management 
inside the stack. On the other hand, the air is introduced to the 
system via an air compressor. Air is pressurized up to the system 
operating pressure after that the temperature of air is adjusted 
to the operating temperature at the heat exchanger in which the 
outlet water stream from the fuel cell is used as coolant. The 
temperature of the fuel cell stack is maintained by a closed cool- 
ing system. Some of the total heat produced in the fuel cell stack 
is assumed to be lost by convection and radiation from the fuel 
cell stack. The remaining heat is taken by the cooling system 
and the outlet streams. 


2.4. Modelling 


2.4.1. Reversible cell voltage 

The reversible cell voltage is the maximum voltage that the 
cell can be produce without any overpotentials and irreversibil- 
ities. There are numerous reversible cell voltage calculation 
equations developed in various forms of the Nernst equation. 
The equation developed by Amphlett et al. [17] is used since 
the equation depends on the data from the Ballard and is a func- 


tion of environment temperature, system pressure and the partial 
pressures of both hydrogen and oxygen. 


Vievy = 1.229 — 8.5 x 1074(Tpc — 298.25) + 4.3085 x 1075 


x Tre [In(pin) + 3 mpo) (9) 


The operating cell voltage is less then the reversible cell 
voltage because of the irreversibilities and overpotentials. The 
hydrogen and oxygen partial pressures may be calculated as: 


1 — xH,0,A 

= P 10 

P = 14+ A/D F Ea/ Ea D) A o 
1 — xH,0,c 

= P, 11 

ree TF A e ai 


where xy,0 is the water mole fraction that is Psat/P, respectively, 
both for anode and cathode. xa and xc are anode and cathode 
dry gas mole fractions, respectively, and £4 and £c are anode 
and cathode stoichiometries, respectively [18]. 


2.4.2. Operating cell voltage 

According to the reversible cell voltage, it is expected that 
the cell voltage remains constant during the operation of the 
fuel cell, but there are some overpotential losses that causes a 
voltage drop with respect to the current density. So the operating 
cell voltage may be expressed by: 


Voperating = Vrev — Virrev (12) 


The overpotentials are mainly classified as: activation losses, 
ohmic losses, and concentration losses. 


Virrev = Vact + Vohm + Veon (13) 
2.5. Activation overpotential 


Activation losses are caused by low reaction rates in both 
anode and cathode by losing some of the energy while driving 
the reactions for transferring electrons. The relationship between 
the current density and overpotential is logarithmic. In this study, 
the activation losses are calculated by using the Tafel equa- 
tions proposed by Barbir and Gomez [9] and Bard and Faulkner 
[19]. The overvoltage at the surface of an electrode follows a 
similar pattern in a great variety of electrochemical reactions 
[20]: 


RTpc i 
Vact, Anode = In | — (14) 
aank io 
RTgc i 
Vact, Cathode = ——> In (=) (15) 
acnF io 


where i is current density (A cm7?); ig is exchange current den- 
sity (A cm~?); R is the universal gas constant (J (kmol K)~}); 
n is the number electrons involved; F is the Faraday’s constant 
(Cmole~!); œa and ac are the empirically determined elec- 
tron transfer coefficient of the reaction at the electrodes at the 
anode and cathode. The exchange current density is higher if 
the reaction rate increases and can be considered as the cur- 
rent density at which the overvoltage begins to move from zero. 
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The ip exchange current density can be determined by empirical 
equation as derived by Berning [21]: 


io(T) = 1.08 x 107?! exp(0.086 x Trc) (16) 


2.5.1. Ohmic overpotential 

Ohmic overpotential arises due to the electrical resistance 
inside the cell. The size of the voltage drop is proportional to the 
current. The resistance to the flow of protons in the membrane 
is much more effective than the resistance of the electrodes, 
resistance of the bipolar plates. The resistance to the flow of 
protons in the membrane is defined as the reciprocal of the proton 
conductivity in the membrane. 


Vohm = 7 Rohm (17) 
t 

Rohm = —— (18) 
Omem 


where Omem is the membrane conductivity (1 Q-!cm7!).The 
membrane conductivity depends on the temperature and also on 
the water content of the membrane as: 


1268 1 
Omem = (0.005139Amem — 0.00326) exp Ee — zzl (19) 
where Amem is the membrane water content. The membrane 


water content is taken from Zawdonowski et al. [22]: 


x J 0043 + 17.8la — 39.85a° + 39.85a°, O<ax<1 
mm | 144+ 14a- 1), 1<a<3 
(20) 
xmoP 
—— (21) 
Psat 


where a is the membrane water activity and xp,o is the water 
mole fraction. 


2.5.2. Concentration overpotential 

Because of the mass transfer limitations and the concentration 
drops at high current densities, some voltage drops occur in both 
anode and cathode. As the reaction proceeds, the concentration 
varies at the surface of the reaction sites and this will create 
a concentration gradient between the reaction sites and bulk 
phase. At high current densities, the concentration gradient will 
be very high and will limit the rate of reactions. Concentration 
overpotential is determined by the equation given by Guzzella 
[23]. 


i Bo 
Uconc = iG ; ) (22) 


Tmax 


where 6; and 62 and imax are constants which depend on the 
operating temperature and reactant concentration. imax and f2 
are taken as 2 also 6, is calculated by the empirical equation of 
Pukrushkapan et al. [24]: 


2.5.3. Fuel cell stack 

Here, in this particular case, the fuel cell stack is composed 
of 97 cells of each having a 900 cm? effective area. The power 
produced by the stack can be calculated by: 


Wstack = Voperating! Acell’cell (24) 


where neel is the number of fuel cells inside the stack; Ace is 
the area of the each cell; i is the current density. The remaining 
exergy calculation for the fuel cell stack is done similar as the 
previous single fuel cell calculations. 


2.5.4. Overall system 
In order to determine the thermodynamic characteristics of 
the system, some general assumptions are made as follows: 


e Both the PEM fuel cell and the fuel cell engine are assumed 

to be at steady state. 

The flows of reactants are incompressible and laminar. 

The gases are assumed as ideal. 

Hydrogen is stored at 10 bar and 298 K. 

Isentropic efficiencies are taken as 70% for compressor, cool- 

ing fan and cooling pump. 

e 20% of the produced heat is assumed to be lost. 

e The relative humidity of the inlet air and hydrogen is taken as 
90%. 

e The base operating parameters are 353 K temperature, 3 atm 
pressure, 0.016cm membrane thickness, 1.5 for both anode 
and cathode stoichiometry. 

e The cost for hydrogen is taken as US$ 10 GW [8]. 

e The capital costs of the system equipments are taken from 
Carlson et al. [25] as US$ 108 kW~! for overall system 


The overall exergy balance is: 


Tec = Ss Ë masin = 5 Easson = 5 Eveat = 5 Ewok 


In order to calculate the system efficiency, the net power 
production is calculated as: 


Wnet = Wec = Weomp ae Weool,pump,act m Wean, act (25) 
Wnet 
Nsys,energy = E (26) 
2 “Ho, in 
Wnet 
Nsys,exergy = e (27) 
in 


where Ejn is calculated from the exergies of the inlet streams of 
the system. 

The overall exergetic cost balance can be written by using 
Eq. (4) as: 


S > Bin Cin) + Ztor = X (Eout,iCout,i) + WnetCw (28) 


we 2 _ 4p. <2, (7.16 x 1074Tpe — 0.622) ( ew Past) + (=1.45 x 10-3 Tae + 1.68) 
0.1173 0.1173 
i= a (23) 
else, (8.66 x 10-5 Tac — 0.068) (z a Pat) + (—1.6 x 10-4 + 0.54) 
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Fig. 2. Polarization curve for fuel cell engine system at different temperatures. 


3. Results and discussion 


In order to investigate the fuel cell system performance, a 
parametric study is conducted. The operating parameters, such 
as temperature, pressure and stoichiometry ratio are considered 
for performance evaluation. 

Fig. 2 shows how the polarization curve of the fuel cell 
changes with current density at various operating temperatures. 
The open-circuit voltage of the system lies between 0.85 and 
0.9 V, primarily depending on the operating temperature. The 
operating parameters other than the investigated one for each 
study, in this case others than temperature, is kept constant in 
the best values, such as 3 atm pressure, 353 K temperature, 1.1 
and 1.5 for fuel and air stoichiometry, respectively, and 0.016cm 
membrane thickness. The effects of the overpotentials are clearly 
seen from this figure that the first region is the activation losses up 
to0.25A cm~2, then the ohmic losses section up to 1.6A cm~2, 
finally the concentration and mass transport losses. Itis also seen 
that the irreversibilities greatly depend on the current density as 
at the low current densities the cell voltage is nearly 0.9 V but at 
high values the voltage decreases up to 0.08 V. Fig. 2 also shows 
that with an increase of the operating temperature the cell voltage 
increases since the reversible cell voltage greatly depends on the 
operating temperature in a direct proportion and the overpoten- 
tials are less affected than the reversible cell voltage by operating 
temperature. The results of the present model are very consistent 
with some experimental studies in the literature [26,27]. Fig. 3 
shows the comparison between the experimental results of Wang 
et al. [27] and the present model. It is seen that the trends of 
the efficiencies are very similar and the values matches greatly, 
more precisely there is a maximum error of 3% while the aver- 
age error is about 1%. These small errors may arise from the use 
of empirical equations used and/or the effect of the deviations 
in operating parameters. 

Fig. 4 shows that in low current densities, the efficiencies 
are increasing up to a peak point at 0.05 A cm~? with 55% in 
energy efficiency and 48% in exergy efficiency. The reason for 
this is that the low current densities, the expenses as the com- 
pressor load, cooling pump load and the fan load are low when 


0.58 O 323K, Experimental Wang et. al. (2003) + 
+ 333K, Experimental Wang et. al. (2003) 
© 343K, Experimental Wang et. al. (2003) 
0.53 A 353K, Experimental Wang et. al. (2003) | 
323 K, Present model 
0.48 — —333 K, Present model — Energy Eff. 
= = = 343 K, Present model —— Exergy Eff. 
ee = 353 K, Present model 
> 0.43 
(S) 
S 
2 0.38 
W 0.33 
0.28 
0.23 
0.18 
0 0.5 4 15 
i [A/cm?] 


Fig. 3. Comparison of the present model with experimental data of Wang et al. 
[27]. 


compared to the values at high current densities [1]. So after a 
critical point with an increase of the power withdrawn in the sys- 
tem the efficiencies start to decrease. When the current density 
reaches high values as 1.6 A cm7? the energy efficiency val- 
ues decreases in a rate higher than the exergy efficiency with 
results of nearly at 2Acm~* current density the efficiencies 
of the system decreased lesser than 5%. On the other hand, 
the effect of the temperature on the system efficiency can be 
seen from this graph. As similar to the polarization curve the 
efficiencies increases with increasing temperature and there is 
a difference of 8% for the high current densities means that 
the operating temperature is a crucial parameter for the fuel 
cells. 

The main results for the study are the exergetic costs of the 
produced electricity from the engine including the expenses 
and all the system capital costs. As it is seen at low cur- 
rent densities, the low operating temperatures are cheaper than 
the high temperatures. But when the current densities higher 
than the 0.8 Acm7?, the situation changes and the price of 
the electricity increases exponentially and the low tempera- 
tures becomes more expensive than high temperatures. This 


Efficiency 


—— Exergy Efficiency 


—-—-—- Energy Efficiency 


0 0.4 0.8 1.2 1.6 2 
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Fig. 4. Effect of temperature on the system performance. 
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Fig. 5. The effect of the temperature on the overall costs and power produced. 


situation arises from the effects of the system components 
exergetic costs as the load of the heat exchanger changes for 
different temperatures and depending on the temperature the 
load of the humidifiers are also changes. Also, for the low 
operating temperatures in high current densities the costs are 
much higher because the produced electricity decreases as it 
is seen in the Fig. 5 and the efficiency decreases with the 
decreasing temperature as consistent with the study of Kazim 
[8]. 

The reversible cell voltage, ohmic overpotentials and con- 
centration overpotentials are essentially pressure dependent, and 
these effects show a kind of nonlinear behaviour, so the effect 
of the pressure can only be investigated from these results. 
Increasing pressure increases the system efficiency as shown 
in Fig. 6. At high current densities the effect of pressure is 
more obvious with a difference of 5% in 1.8 A cm~? current 
density but at low values gap between the 2.5 and 4atm can 
be constant for up to 0.9 A cm~? with an efficiency of 3%, 
respectively. 

Although, high pressure operation requires pressurization of 
inlet streams that increases the compressor load, the net power 
production increases with the increasing pressure as it is seen 
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Fig. 6. Effect of pressure on the system performance in terms of energy and 
exergy efficiencies. 
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Fig. 7. Variations of production cost and network output with the current density 
at various operating pressures. 


from Fig. 7. Since the decrease in the overpotentials are higher 
than the increase in the compressor load that brings an increase 
in the overall efficiency with an increasing pressure. As it is 
expected, the effect of the pressure on the cost of the unit elec- 
tricity production is in a positive way as shown in Fig. 7. Since 
the efficiency of the system increases depending on the pres- 
sure, the production cost for the 2.5 atm case is higher than the 
cost for the 4atm case, there is a maximum of US$ 7 GW7! 
difference between the two cases. The results for the effect of 
pressure are also consistent with experimental studies of Wang 
et al. [27]. 

The effect of the membrane thickness for the system per- 
formance is given in Fig. 8. It is seen that with increasing the 
membrane thickness the efficiency of the system decreases down 
to 2% for moderate current densities. As the thickness of the 
membrane increases, the network output also decreases as seen 
in Fig. 9. The effect is much clearer at the high current densi- 
ties with a difference of 4kW. This decrease is expected, since 
the ohmic overpotential losses highly depends on the membrane 
thickness. Of course, it must be considered that the fuel crossover 
effect which decreases the system efficiency at low membrane 
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Fig. 8. Variation of energy and exergy efficiencies with the current density at 
various membrane thicknesses. 
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Fig. 9. Effect of membrane thickness on the production cost and power output. 


thickness is neglected for this study and the fuel utilization 
ratios are kept constant depending on the fuel stoichiome- 
try. Beside the efficiencies and power, the cost for electricity 
production is also affected directly by an increase of the mem- 
brane thickness (Fig. 9). The maximum difference appears to be 
US$ 3 GW- ! for the highest current density values at different 
thicknesses. 

The effect of the cathode stoichiometry, which is the ratio of 
the fed oxygen to the required oxygen that gives the magnitude 
of the excess oxygen, on the system performance, is clearly 
shown in Fig. 10. By the increase of excess oxygen, the effi- 
ciency drops by about 1% since the exergy inlet to the system 
by the oxygen fed is increases but this increase does not causes 
a performance increment for the system so the efficiency of the 
system decreases dramatically. Fig. 11 shows the effect of the 
anode stoichiometry that means the ratio of the fed hydrogen 
to the required hydrogen to the system. Since the hydrogen is 
the main energy (fuel) source of the fuel cell, as expected the 
effect of the anode stoichiometry is very high on the efficiency 
values. With an increase of anode stoichiometry values from 1.1 
to 3, the efficiency drops down to 35% with increasing current 
density. 
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Fig. 10. Effect of cathode stoichiometry on exergy efficiency. 
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Fig. 11. Effect of anode stoichiometry on system performance in terms of energy 
and exergy efficiencies. 


Fig. 12 shows the exergetic cost flows of the system equip- 
ment versus the current density. It is seen that the cost flows 
of the equipment do not depend on the current density. When 
these results are compared to the ones in Fig. 13, as the cost 
flows of the fuel cell stack in different temperatures with respect 
to the current density, it is seen that the cost of stack is much 
higher than costs for other components. These results are in 
agreement with the experimental data obtained by Carlson et 
al. [25]. Such a large difference means that the fuel cell stack 
is the most costly component based on the overall electric- 
ity production cost with having more than 40% contribution 
besides the stack cost is highly dependent to the current den- 
sity since the maximum irreversibilities occur in fuel cell stack 
and these irreversibilities are highly dependent to the current 
density as activation and ohmic overpotentials. Also, the effect 
of the temperature to the fuel cell stack cost can be seen 
from Fig. 13, which is also compatible with the overall system 
cost since increasing temperature decreases the fuel cell stack 
cost. 
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Fig. 12. System components cost vs. current density. 
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Fig. 13. Fuel cell stack cost vs. current density. 


4. Conclusions 


An exergoeconomic model has been developed and applied 
to a PEM fuel cell engine system for transportation applica- 
tions. A parametric study is conducted to investigate the system 
performance and the cost behaviour depending on the operat- 
ing parameters. It is found that, with an increase of operating 
temperature the system efficiency increases and the overall pro- 
duction costs decrease. Also, high pressure is another positive 
parameter for the system efficiency if we neglect the increase 
of manufacturing costs at higher pressures. The increase in the 
cathode stoichiometry leads to a small decrease in the over- 
all system exergy efficiency, but the anode stoichiometry has 
a major effect on the efficiency. The main contribution to the 
overall cost is made by the fuel cell stack that has the high- 
est irreversibility compared to the other system components. So 
any increase in the stack efficiency will greatly affect the over- 
all performance and the production cost, and will contribute the 
commercialization of the fuel cell systems in the transportation 
sector. 
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